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ABSTRACT 


Under  typical  conditions  in  the  subtropical  thermocline 
double-diffusion  occurs  in  the  form  of  salt  fingering 
(warm,  salty  fluid  overlies  cold,  fresh  fluid) .  The 
formation  of  staircases  in  the  thermohaline  structure  of 
the  ocean  has  been  observed  since  the  late  1960s,  with 
recent  field  data  collected  within  the  tropical  Atlantic 
displaying  staircases  with  high-gradient  interfaces 
characterized  by  a  unique  spatial  orientation  determined  by 
background  temperature  and  salinity.  Competing  theories 
have  been  proposed  to  explain  the  dynamics  of  these 
staircases;  however,  the  origin  of  the  staircases  and  the 
mechanism  that  controls  final  equilibrium  remains  poorly 
understood.  This  thesis  examines  staircase  development  in 
the  tropical  Atlantic.  Incorporating  double  diffusion  using 
the  flux-gradient  formulation  of  Radko  and  Smith, 
staircases  are  numerically  simulated  to  resolve  the 
controversial  aspects  of  the  staircase  theories.  The 
staircase  simulations  are  critically  evaluated  against  the 
C-SALT  experiment  observations,  and  the  conditions  for 
their  formation  and  explanation  of  the  fully  equilibrated 
state  are  resolved  and  explained.  The  effect  of  acoustic 
propagation  through  the  three-dimensional  modeled 
staircases  is  evaluated  to  determine  the  impact  of  these 
large  staircase  areas  on  various  frequencies  and  depths. 
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INTRODUCTION 


I . 

A.  DOUBLE  DIFFUSION 

1.  History 

Double  diffusion  is  a  phenomenon  that  occurs  due  to 
the  differing  molecular  dif fusivities  of  two  properties 
(e.g.;  temperature  and  salinity),  within  a  stably 
stratified  fluid.  Double-diffusion  operates  in 

configurations  in  which  these  properties  have  an  opposing 
effect  on  the  stability  of  the  water  column. 

Double  diffusion  is  a  relatively  recent  discovery  in 
the  realm  of  science  (Stern  1960),  although  W.  Stanley 
Jevons,  an  English  economist,  came  very  close  to 
discovering  it  in  1857  when  he  created  sugar  fingers  during 
an  experiment.  Although  Jevons  understood  that  heat 
diffused  faster  than  the  dissolved  sugar,  his  statement 
that  "all  convection  should  occur  as  minute  streamlets" 
failed  to  capture  the  physics  behind  the  phenomenon 
(Schmitt  1995)  . 

Whilst  others  would  go  on  to  attempt  to  explain  the 

phenomena,  including  the  English  physicist  Lord  Rayleigh, 

the  subject  was  correctly  explained  only  in  the  mid-20th 

century  by  Melvin  Stern  (1960) .  Stern's  discovery  was 

motivated  by  an  experiment  by  his  colleagues,  Stommel, 

Arons  and  Blanchard  (1956),  who  developed  the  Perpetual 

Salt  Fountain.  This  was  an  experiment  initially  designed  to 

test  pressure  at  the  base  of  an  ocean,  by  inserting  a 

narrow  heat-conducting  pipe  vertically  to  draw  water  to  the 

surface.  Due  to  the  lower  salinity  of  the  basal  water,  once 

at  thermal  equilibrium  it  would  continue  to  rise  towards 

1 


the  surface  due  to  being  lighter  than  its  surrounds.  While 
investigating  this  idea,  the  concept  that  placing  a  spigot 
on  the  top  of  the  tube  would  cause  the  water  to  run 
continuously  was  introduced,  thus  the  "Perpetual  Salt 
Fountain"  was  born.  Stern  took  this  experiment  to  the  next 
step,  realizing  that  salt  fountains  occur  naturally  in  the 
ocean,  as  a  result  of  heat  conducting  100  times  faster  than 
the  diffusion  of  salt  in  the  ocean.  Hence  thermal 
eguilibrium  of  a  parcel  of  water  with  its  surroundings  can 
occur  faster  than  chemical  equilibrium,  resulting  in  an 
ocean  containing  parcels  with  differing  salinities  adjacent 
to  each  other  devoid  of  a  salt  fountain  pipe.  Different 
salinities  result  in  different  density  gradients,  and  thus 
motion  dubbed  as  double  diffusive  convection. 

2 .  Salt  Fingers 

There  are  two  types  of  double  diffusive  convection  in 
the  ocean.  The  first  is  when  cold,  fresh  water  overlies 
salty,  warm  water.  If  a  warm  salty  parcel  is  displaced 
upward,  due  to  the  difference  in  molecular  diffusivity  it 
cools  much  quicker  then  it  loses  salinity,  thus  becoming 
denser  than  its  surroundings  and  sinking  back  down.  As  it 
sinks  it  regains  heat  quickly,  but  is  now  less  saline  than 
its  surroundings,  thus  rising  again,  this  time  further  than 
it  did  previously.  This  process  can  continue  to  occur, 
slowly  mixing  the  water  column.  This  is  known  as  diffusive 
convection.  The  other  type  is  salt  fingering,  which  occurs 
when  warm,  salty  water  overlies  cold,  fresh  water.  If  the 
warm,  salty  parcel  is  displaced  downward,  as  in  the 
previous  example  it  will  lose  heat  quicker  than  salinity, 
becoming  denser  than  its  surroundings,  and  thus  continuing 
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to  sink.  These  parcels  appear  as  thin  elongated  streams, 
only  1  to  5  centimeters  wide,  and  tens  of  centimeters  long. 


Salt  Fingers 
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Figure  1.  Schematic  diagram  illustrating  the  physical 
mechanisms  of  salt  fingering  (From  Ruddick  and 

Kerr,  2003) 

3 .  Thermohaline  Staircases 

In  the  ocean,  vertical  stratification  is  often 
characterized  by  the  appearance  of  well-mixed  layers 
separated  by  fine  interfaces  with  large  temperature  and 
salinity  gradients,  known  as  thermohaline  staircases.  The 
interfaces  of  these  staircases  are  occupied  by  salt 
fingers,  which  cause  a  net  unstable  buoyancy  flux  that 
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stirs  the  layers  (Stern  and  Turner  1969)  .  This  stirring  of 
the  layers  can  result  in  vertical  growth  of  the  layer 
through  destruction  of  the  interface  and  merging  of  the 
layers . 

There  is  a  large  abundance  of  these  staircase  regions 
throughout  the  world,  and  the  increased  temperature  and 
salinity  fluxes  in  these  staircases  (Schmitt  2005)  suggest 
that  they  play  an  important  role  in  global  ocean  mixing. 


4 .  Theoretical  Background 


The  density  ratio  controls  the  intensity  of  salt 
fingering  and  is  defined  by: 


_  c/AT 
p  PAS 


(1) 


where,  a  is  the  thermal  expansion  coefficient  and  (3  is  the 
saline  expansion  coefficient.  If  the  density  ratio  is  less 
than  one,  then  the  density  gradient  is  unstable  and  will 
overturn.  High  values  of  density  ratio  correspond  to  weak 
double  diffusive  mixing,  whereas  density  ratios  that  are 
slightly  above  unity  are  associated  with  vigorous 
fingering . 


Another  important  parameter  in  the  formation  of  salt 
fingers  is  the  flux  ratio. 


7  = 


aFT 

Wn 


(2) 


where  FT  and  Fs  are  the  molecular  fluxes  of  heat  and  salt, 
respectively.  This  value  must  theoretically  be  below  one  to 
satisfy  the  net  transport  of  density  "up  gradient." 
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B.  C-SALT 

In  1985,  a  field  program  by  the  name  of  Caribbean- 
Sheets  and  Layer  Transects  (C-SALT)  was  carried  out  in  the 
tropical  North  Atlantic  Ocean  east  of  Barbados  to  survey 
the  thermohaline  staircases  in  the  area,  developed  due  to 
vertical  mixing  by  salt  fingers.  The  survey  area  lies 
between  150m  and  750m  depth,  at  the  confluence  of  the 

highly  saline  Subtropical  Underwater  and  the  fresher 
Antarctic  Intermediate  water,  and  this  gives  rise  to  a  very 
strong  destabilizing  salinity  gradient  which  can  drive  salt 
fingers  (Schmitt  1987)  .  The  temperature  profiles  in  the 

survey  were  classified  as  layered,  weakly  layered  or 

nonlayered,  with  an  area  within  10°  and  14°N,  and  52°  and 
57°W  having  the  thickest  mixed  layers,  largest  number  of 
layers,  and  largest  temperature  differences,  as  seen  in 

Figure  2 . 


5 


II.  LATERAL  COHERENCE,  MERGING,  AND  ALIGNMENT  OF 
THERMOHALINE  STAIRCASES 


A .  LATERAL  COHERENCE 

1.  C-SALT  Observations 

An  interesting  feature  of  the  C-Salt  observations  was 
lateral  coherence  of  the  staircase  layers  over  many 
kilometers  by  examining  profiles  at  a  variety  of 
separations.  The  shortest  separation  between  consecutive 
profiles  was  1.3  km  using  a  CTD  tow-yo,  covering  a  total 
distance  of  33  km.  The  results  of  this  are  shown  in  Figure 
3,  whereby  individual  layers  can  be  traced  over  the  entire 
hydrographic  section. 


TEMPERATURE 


t»  9  12  15 


TEMPERATURE 

Figure  3.  Three-dimensional  view  of  CTD  tow-yo  (rear 

profile  to  the  south)  as  the  ship  steamed  to  the 
north.  1.3  km  separation  between  profiles,  33  km 
total  distance  covered  (From  Schmitt  et  al .  1987) 
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This  lateral  coherence  was  not  limited  to  33  km,  as  it 
was  also  measured  in  individual  layers  over  a  distance  of 
more  than  100  km  using  aerial  dropped  bathythermographs 
(XBT)  spaced  every  5.5  km  in  an  east  to  west  pattern 
(Figure  4 )  . 


TEMPERATURE  (  Or*  t  ) 
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Figure  4.  XBT  temperature  profiles  from  east  to  west,  with 
5.5  km  separation  between  drops  and  total  distance 
of  220  km.  The  profile  at  far  left  is  correct,  each 
subsequent  profile  is  offset  by  1.6  °c  (From 
Schmitt  1987) . 

Over  longer  distances  there  were  significant 
temperature  and  salinity  changes  within  the  layers,  however 
an  examination  of  the  temperature  and  salinity 
characteristics  of  these  layers  produced  a  surprising 
pattern  where  temperature  and  salinity  values  are  grouped 
in  distinct  lines,  suggesting  layers  get  warmer,  saltier 
and  denser  to  the  north  and  west  (Schmitt  1988)  .  Figure  5 
displays  this  unique  grouping,  which  demonstrates  that  the 
well-defined  slope  of  the  layers  crosses  isopycnals.  The 
horizontal  density  ratio,  when  averaged  over  the  ten 
principle  layers,  was  0.85  (±0.02)  in  the  spring,  and  0.84 

(±0.03)  in  the  fall.  Among  other  things,  this  value  aided 
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in  the  realization  that  salt  fingers  were  the  dominant 
process  in  maintaining  C-SALT  staircases,  as  the  value  of 
the  lateral  density  ratio  was  distinctly  different  from 
vertical  mixing  by  small-scale  turbulence  (R  —1.6),  and 
isopycnal  eddy-stirring  (1.0)  (Schmitt  1988). 


"*  -*-* - r*~  ■»  —  •  |  — i —  •  - l—  - 1 — 
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SALINITY  (  PSU  ) 

Figure  5.  Potential  temperature— salinity  values  of  layers 
from  CTD  stations  during  C-SALT.  Solid  circles  are 
from  layers  greater  than  10  m  thick,  open  circles 
are  from  layers  5-10  m  thick.  Hashed  lines  are 
temperature-salinity  relationship  obtained  from 
stations  at  northwest  and  southeast  corners  of  the 
survey.  Solid  lines  represent  isopycnal  surfaces 
(From  Schmitt  1987) . 
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B .  MERGING 

1 .  Theory 

It  is  difficult  to  gauge  the  evolution  of  layers  in 
the  C-SALT  observations,  short  of  saying  that  some 
staircases  have  persisted  over  15  years,  having  again  been 
observed  in  2001  (Schmitt  2005)  .  The  evolution  and  long 
term  persistence  of  these  staircases  can  only  be 
speculated,  and  the  question  remains  as  to  whether  they 
form  over  decades  or  centuries,  and  what,  if  anything, 
causes  a  collapse  in  the  staircase.  The  theory  behind 
staircase  layer  merging  is  explored  here,  and  validated 
using  a  numerical  model  discussed  later. 

It  has  been  proposed,  and  confirmed  by  direct 
numerical  simulation  (DNS) ,  that  layer  formation  occurs  as 
a  result  of  the  instability  of  the  heat/salt  flux  ratio,  y, 
that  occurs  when  the  flux  ratio  is  a  decreasing  function  of 
density  ratio  (Radko  2003) .  Further  to  this  proposal  is  the 
theory  that  relatively  thin  layers  within  the  staircase 
merge  until  the  thickness  of  the  layers  exceeds  a  critical 
value,  H0,  with  equilibrium  occurring  as  a  result  of 
inhomogeneity  of  the  convecting  layers  stabilizing  the 
staircase  (Radko  2005) . 

In  Radko  (2007),  a  merging  theorem  was  developed  which 
states : 

If  the  buoyancy  flux  decreases  with  the  buoyancy 
variation  across  the  step,  then  the  staircase  is  affected 
by  the  B-instability,  and  the  spontaneous  merger  of  layers 
occurs  through  gradual  erosion  of  weak  interfaces.  If, 
however,  fluxes  increase  with  step  height,  then  coarsening 
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of  a  staircase  may  result  from  drift  and  collision  of  the 
adjacent  interfaces.  [Radko,  2007]. 

In  Chapter  III,  we  analyze  the  merging  events  in  the 
numerical  model,  using  the  initial  temperature  of  a  layer 
at  commencement  of  a  merging  event  and  the  temperature 
variation  in  the  layer  during  the  event.  In  all  simulations 
analyzed  in  this  study,  mergers  occur  following  the  13- 
merger  scenario. 

C.  ALIGNMENT  WITH  BACKGROUND  FLOW 

Another  area  of  research  extending  from  the  C-SALT 
experiment  that  was  explored  further  by  Wall  (2007)  was  the 
alignment  of  the  slope  of  staircase  interfaces  with  the 
background  large  scale  isotherms  within  the  staircase 
region.  Wall  first  began  by  considering  a  mixed  layer 
sandwiched  between  two  interfaces  as  shown  in  Figure  6,  and 
starting  with  the  conservation  law  for  temperature  and 

salinity  in  a  mixed  layer  ( FT2  and  FTl  are  the  temperature 
flux  at  the  top  and  bottom  of  the  layer,  respectively) . 
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Figure  6.  Basic  two  layer  interface  (from  Wall  2007) 

By  adopting  a  co-ordinate  system  that  is  in  the  same 
direction  as  the  flow,  and  using  the  Radko  (2005) 
hypothesis  that  as  a  staircase  evolves  in  time  to  a  quasi¬ 
equilibrium  state,  the  interfacial  flux  ratio  approaches 
its  minimum  value,  ymin  ,  Wall  arrived  at  the  following: 

(jjSa  ' 3x)  =  -  Fo  =  (3) 

0s12  /  3x)  (fs2  -  fsi) 

Wall  rearranged  the  left  hand  side  of  this  equation  to  be 

c)T 

^l)S^L  anc*'  determined  this  value  from  the  temperature- 

salinity  plot  similar  to  that  shown  in  Figure  5  for 
stations  in  the  C-SALT  area  exhibiting  strong  staircases. 
The  temperature-salinity  variation  along  the  mixed  layer 
was  derived  in  terms  of  Cartesian  derivatives  and  after  a 
series  of  derivations  Wall  produced: 
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(4) 


.  dT 

tga  =  (~tgaT  —  +  yn 
dz 


dS 

*tgas  —)/(7n 
dz 


tga  =  (- tgaT  *  Rp  +  7min  *  #as)/  (7,, 


dS_dT) 

dz  dz 

n-Rp) 


(5! 


where  aT  and  as  are  the  slope  of  the  background  temperature 

and  salinity,  respectively,  a  is  the  slope  of  the  layer 
interface,  and  assuming  Rp>Ymin  yields: 

tga  ~  tgaT  ( 6 ) 

Wall  thus  hypothesized  that  the  interface  should  align  with 
isotherms  but  not  necessarily  isohalines.  An  observational 
test  was  carried  out  on  the  C-SALT  staircases  discussed 
above,  using  a  numerical  model  based  on  the  GFDL  M0M3  code 
(Pacanowski  and  Griffies  1999)  to  determine  the  direction 
of  the  mean  background  flow.  The  results  of  this 
observational  test  showed  a  weak  relationship  in  the 
direction  of  background  flow  between  the  isotherm  and  the 
slope  of  the  layer  (much  less  so  in  the  case  of 
isohalines) .  In  the  direction  normal  to  background  flow, 
the  relationship  between  the  isotherm  and  slope  of  the 
layer  disappeared.  This  idea  is  further  explored  in  the 
following  chapter,  in  order  to  attain  a  stronger 
confirmation  of  this  feature  shown  by  Wall. 
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Ill .  SIMULATION  OF  STAIRCASES-EVOLUTION  AND 

ALIGNMENT 


A .  MODEL  SETUP 

1.  Massachusetts  Institute  of  Technology  General 
Circulation  Model  (MITgcm) 

In  order  to  simulate  the  staircases  seen  during  the  C- 
SALT  field  program,  the  MITgcm  numerical  model  was  chosen. 
The  desirable  features  of  this  coupled  model  include,  but 
are  not  limited  to,  the  non-hydrostatic  capability  that 
allows  the  model  to  simulate  fine,  meso,  and  synoptic  scale 
processes,  and  the  ability  to  perform  on  a  wide  variety  of 
computational  platforms  (Adcroft  et  al.  2012)  .  The 
flexibility  in  building  the  model  with  varying  parameters 
was  ideal  for  the  purpose  of  this  investigation,  where  many 
parameters  such  as  the  horizontal  and  vertical  diffusivity 
rates  are  not  known,  but  can  be  slowly  adjusted  to  produce 
simulations  that  reflect  conditions  in  nature. 

Model  runs  were  performed  on  the  Department  of  Defense 
High  Performance  Computing  Modernization  Program,  Cray  XE6 
Open  Research  System  (Chugach) ,  located  at  ERDC,  Vicksburg, 
Mississippi.  128  processors  were  used  in  coarse  resolution 
grid  simulations  encompassing  the  entire  ocean  basin. 

During  early  model  initialization  and  simulation  a 
coarse  uniform  1  degree  grid  resolution  was  used  over  the 
entire  ocean  basin.  Once  this  model  had  reached  a  period  of 
perceived  equilibrium  at  time  =  18  yrs,  a  high  resolution 
(Ax  =  Ay  =  10  km)  interpolation  was  carried  out  in  the  area 
of  staircases,  with  the  remainder  of  the  ocean  basin 
staying  at  coarse  grid  resolution.  The  simulation  used  Az=l 
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m  and  a  Cartesian  grid.  Initial  runs  involved  the  use  of 
the  Gent  McWilliams /Redi  SGS  eddy  parameterization  scheme 
(GMredi)  which  mixes  tracer  properties  along  isentropes 
(Redi  1982),  and  adiabatically  rearranges  tracers  through 
an  advective  flux  where  the  advecting  flow  is  a  function  of 
the  slope  of  the  isentropic  surface  (Gent  and  McWilliams 
1990) .  The  GMRedi  scheme  was  implemented  using  the  tapering 
scheme  of  Large  et  al .  (1997);  however,  it  was  found  that 
over  time  the  vertical  temperature  salinity  profiles  were 
smoothed  to  such  an  extent  that  the  staircases  were 
removed,  an  unexplained  result  in  the  simulation  in  direct 
opposition  to  what  we  know  happens  in  the  ocean.  Thus  the 
GMredi  scheme  was  not  used  in  most  simulations. 

Trial  simulations  were  also  carried  out  using  various 
values  of  biharmonic  diffusion  coefficients,  in  order  to 
accurately  account  for  both  mixing  as  a  result  of 
turbulence  and  double  diffusion.  The  biharmonic  diffusion 

coefficient  was  ultimately  parameterized  at  1  * 1 0  6  m4 1  s  over 
all  1000  vertical  layers. 

Various  parameter  adjustments  were  made  throughout  the 
course  of  model  experimentation,  and  these  are  discussed 
later  in  the  chapter. 

2 .  Initial  Datasets 

At  model  initialization,  temperature  and  salinity 
datasets  were  generated  to  provide  a  density  ratio  pattern 
over  the  tropical  north  Atlantic  qualitatively  resembling 
observations  during  C-SALT.  The  initial  modeling  of 
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temperature  and  salinity  used  in  this  experimentation  at 
the  surface  and  1000m  is  shown  along  with  density  ratio 
values  in  Figure  7 . 
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Figure  7.  Temperature  and  Salinity  at  the 
1000m  depth,  and  Density  Ratio,  at 
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surface  and  at 
initialization . 


3. 


Model  runs  performed  without  Double  Diffusion  and 
with  Constant  Flux  Ratio 


In  order  to  gauge  the  effect  of  double  diffusion  on 
the  simulation,  several  model  runs  were  conducted  using 
different  K  profile  parameterization  (KPP)  setups.  The  KPP 
scheme  within  MITgcm  is  a  product  of  Large  et  al .  (1994), 

and  it  unifies  the  treatment  of  a  variety  of  unresolved 
processes  in  vertical  mixing  (Adcroft  et  al .  2012)  .  The 

functionality  to  switch  double  diffusion  on  and  off  within 
the  KPP  scheme  allows  the  ability  to  display  a  profile 
within  the  expected  staircase  region  with  and  without 
double  diffusive  mixing,  so  as  to  see  the  effect  double 
diffusion  (and  in  this  case,  salt  fingers)  has  on  the 
development  of  staircases.  Another  useful  feature  of  the 
KPP  scheme  in  MITgcm  is  the  ability  to  adjust  the  flux 
ratio.  The  flux  ratio  model  used  for  all  modeling  is 
suggested  by  Radko  and  Smith  (2012),  who  used  high 
resolution  DNS  data  to  obtain  the  following  algebraic 
expression : 


Y  =  a„  exp  (-b  R  )  +  c 

1  Y  “  '  y  p'  Y  ( 7  ) 

where  (a^,  by,  cy)  =  (2.709,  2.513,  0.5128).  To  illustrate  the 

effects  of  using  a  variable  flux  ratio,  another  simulation 
was  also  conducted  with  a  constant  flux  ratio  of  0.5128.  It 
is  suggested  in  Radko  and  Smith  that  for  the  step-like 
patterns  of  thermohaline  staircases  to  occur,  the  flux 
ratio  must  decrease  with  an  increasing  density  ratio,  in 
the  oceanographically  relevant  range  (1<  Rp<3)  .  We  put  this 
theory  to  the  test  by  comparing  these  two  simulations  with 
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the  expectation  that  with  a  constant  flux  profile  a 
smoothed  vertical  temperature  profile  would  be  evident,  as 
opposed  to  a  flux  ratio  utilizing  Equation  10,  where  a 
stepped  vertical  temperature  profile  would  be  expected. 
Figures  8  through  10  show  a  comparison  of  the  varying  KPP 
examples  discussed  here  at  time,  t=18  years,  and  it  can  be 
seen  that  the  ideal  KPP  set  up  for  growth  and  evolution  of 
staircases  similar  to  those  observed  in  C-SALT  is  using  the 
flux  ratio  at  Equation  7  with  double  diffusion  as  a 
dominant  vertical  mixing  process.  Profiles  are  taken  at  the 


location  x  =  —Lr  and  y  =  —L 

3  4  y 


but  the  patterns  are 


representative  of  the  stratification  in  the  entire  domain. 


Vertical  Temperature  profile  with  constant  flux 
ratio,  y  =  0.5128  . 
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Figure  8 . 


Depth  (m)  hi  Depth  (m) 


9.  Vertical  Temperature  profile  with  the  effects  of 
double  diffusion  removed  from  the  model. 


Figure  10.  Vertical  Temperature  profile  with  flux  ratio  from 

Equation  10  and  active  double  diffusion. 
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B. 


MODEL  RESULTS 


1.  Effects  of  Vertical  Eddy  Dif fusivities 

Vertical  eddy  dif fusivities  due  to  small-scale 
turbulence,  KT  and  Ks  ,  could  be  specified  within  MITgcm 

and  typical  oceanic  values  are  of  the  order  1CTJ  m2  /  s  .  In 
order  to  examine  the  impact  of  turbulent  dif fusivities  on 
thermohaline  layering,  the  model  was  assigned  a  series  of 
KT  and  Ks  values  ranging  from  0  to  10-4  m2/s  ,  and  then  run 
to  t~18  years  on  each  simulation.  A  Matlab  code  was 
developed  to  compare  the  staircase  distributions  with  a  two 

dimensional  map  of  the  density  ratio,  R  ,  between  750  and 

850  m.  To  identify  staircase  distributions,  the  average 
density  ratio  over  the  100  m  vertical  section  was 
calculated,  and  areas  where  density  ratio  changes  within 
this  area  were  greater  than  10%  of  the  average  were 
identified: 


R 


p 


max  R  -  min  R 

- 5 - 2-  >  AR 

1000  p 


(8) 


where  Rp  is  10%  of  the  average  density  ratio  and  ARp  is  the 

change  in  density  ratio.  If  the  individual  density  ratio 
change  between  750  and  850  m  was  greater  than  10%  of  the 
average  density  ratio  change,  and  more  than  a  third  of  the 
vertical  profile  (between  750  and  850  m)  at  individual 
points  exhibited  this  characteristic,  than  staircases  were 
deemed  to  exist  at  this  location. 


The  results  are  displayed  in  Figures  11  through  14, 
where  the  white  contours  indicate  the  presence  of 
thermohaline  staircases.  Alongside  the  staircase  area  is  a 
graphical  representation  of  the  flux  ratio  dependence  on 
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the  density  ratio  at  the  specified  eddy  diffusivity.  As  per 
the  suggestion  by  Radko  and  Smith  (2012),  thermohaline 
staircases  are  prevalent  when  the  flux  ratio  decreases,  but 
are  almost  non-existent  when  the  flux  ratio  increases  with 
increasing  density  ratio.  It  is  evident  that  when  KT= 2.5* 

1CT6  m2  /  s  and  ^S=2.5*1CT6  m2  /  s ,  the  numerical  model  closely 
replicates  the  observed  density  ratios  and  staircase 
formation  areas  of  the  C-SALT  experiment,  and  at  these 
regions  the  flux  ratio  is  also  decreasing.  Other  parameters 
that  were  adjusted  to  best  replicate  C-SALT  observations 
are  in  Table  1 . 


Parameter 

Value 

Vertical  Eddy  Viscosity 

mo-4  m2/s 

Lateral  Eddy  Viscosity 

2*  10 5  nils 

Lateral  Eddy  Diffusivity 

1*102  nils 

Table  1.  Parameter  values  in  the  numerical  model. 
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R(  Longitude  (deg) 

Figure  11.  Total  flux  ratio  versus  density  ratio  (left)  and 
density  ratio  (right)  for  Kt=Ts=  0  m1  / s  .  Upon 
analysis  staircases  are  evident  throughout  the 

region . 


Figure  12.  Total  flux  ratio  versus  density  ratio  (left)  and 
density  ratio  (right)  for  KT=KS  =  1CT6  mis. 
Staircases  are  largely  prevalent  throughout  the 
region  as  indicated  by  contours. 
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Figure  13.  Total  flux  ratio  versus  density  ratio  (left)  and 
density  ratio  (right)  for  KT=KS=  2.5*10^  mis. 
Staircases  are  evident  a  region  similar  to  that 
observed  in  the  C-SALT  experiment. 


Figure  14.  Total  flux  ratio  versus  density  ratio  (left)  and 
density  ratio  (right)  for  KT=KS  =  105  mis. 
Staircases  are  almost  completely  absent  at  this 
vertical  eddy  diffusivity. 
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2 .  Mechanism  of  Merging  and  Growth  Rates 

The  model  output  was  examined  from  initiation  to 
perceived  equilibrium.  The  images  in  Figure  15  through  17 
represent  a  three  dimensional  visualization  of  the 

staircase  looking  north  east  from  50  West  and  8  North 
between  depths  of  500  to  800  m  shown  at  initialization  of 
the  model,  an  early  point  in  the  merging  stage  of  the 
staircase,  and  fully  matured  merged  layers.  This  three 
dimensional  rendering  is  produced  by  VAPOR 

(www.vapor.ucar.edu),  a  product  of  the  Computational 
Information  Systems  Laboratory  at  the  National  Center  for 
Atmospheric  Research. 


Figure  15.  Three-dimensional  visualization  at  t=0  secs. 
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Figure  16.  Three-dimensional  visualization  at  t=10  years. 

Staircases  are  evident  at  this  time  and  layers  have 
commenced  growth  through  merging. 


Figure  17.  Three-dimensional  visualization  at  t=18.6  years 
with  evident  staircases  merged  to  equilibrium. 
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The  evolution  of  the  temperature  profile  in  time  is 
shown  in  Figures  18  and  19,  from  the  initialization  (t=0)  , 
the  time  t=5 . 4  years  at  which  flux  ratio  instability 
modifies  the  initial  profile,  the  time  t=9.2  years  where 
numerous  density  inversions  are  evidenced  within  the 
profile,  to  the  point  where  a  well-defined  staircase  is 
evident  at  t=10.4  years,  and  then  several  time  steps 
illustrating  the  merging  of  the  steps  through  to  perceived 
equilibrium  at  t=18.6  years.  Profiles  are  taken  at  the 


location  x  =  —Lr  and  y  =  —L,  or  50 
3  '  4  -v 


W,  8  N  on  a  global 


grid , 


Figure  18.  Formation  and  evolution  of  layers  in  the 

numerical  model  shown  at  (a-d)  t=0,  5.4  ,  9.2,  10.4 

years . 
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Figure  19.  Evolution  of  layers  in  the  numerical  model  shown 

at  (a-d)  t=14.7,  15.9,  17.7,  18.6  years. 

The  temperature  gradient  as  a  function  of  time  is 
displayed  in  Figure  20,  revealing  the  merging  events  over 
time.  From  this  image  it  is  evident  that  as  the  strong 
interfaces  grow  the  weaker  interfaces  gradually  decay  and 
disappear  (strong  interface  growth  at  the  expense  of  the 
weaker  interface,  the  B-merger  as  discussed  in  Chapter  II) . 
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Figure  20.  Space-time  diagram  of  the  numerically  modeled 

temperature  gradient  (  elm  ) .  Layers  merge  when  the 
relatively  strong  interfaces  grow  further  at  the 
expense  of  weaker  interfaces,  which  decay  and 
eventually  disappear  (B-merger  scenario) . 


Essential  insights  into  the  layer-merging  process 
brought  by  theoretical  studies  (Huppert  1971  and  Radko 
2005)  ,  is  that  the  merging  effect  is  a  form  of  linear 
instability.  Thus  perturbations  for  a  quasi-steady  state 
are  expected  to  grow  exponentially.  Our  simulations  offer  a 
convenient  opportunity  to  test  these  theories.  Figures  21 
and  22  show  two  separate  merging  events  selected  from  the 
DNS  over  the  time  frame  of  the  merge.  The  growth  rate  is 
exponential  for  both  examples,  which  supports  the  theory 
merging  as  a  result  of  linear  instabilities. 
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Figure  21.  Change  in  temperature  vs.  time  during  merging 

event  1 . 


Figure  22.  Change  in  temperature  vs.  time  during  merging 

event  2 . 

30 


3. 


Lateral  Coherence  of  Modeled  Layers 


Another  feature  of  C-SALT  that  was  also  evidenced  in 
the  MITgcm  simulation  was  lateral  coherence  of  the  layers 
over  many  hundreds  of  kilometers.  The  lateral  extent  of 
these  layers  within  the  simulation  measured  greater  than 
300  km  with  minor  spatial  changes  in  temperature  and 
salinity  within  the  layers  as  was  also  evidenced  in  C-SALT. 
Temperature  and  salinity  data  from  three  individual  layers 
is  plotted  in  Figure  23,  and  the  result  clearly  portrays 
the  layer  temperature  and  salinity  grouping  aligned  in  a 
straight  line  with  the  layer  slope  crossing  isopycnals,  as 
it  did  during  C-SALT  (Figure  5)  .  It  is  noted  that  this 
pattern  became  more  pronounced  as  horizontal  eddy 
dif fusivities ,  KhT  and  KhS  decreased  by  0(10),  from  1000  to 

100  m1  / s  .  The  average  density  ratio  was  defined  by  the 
temperature  and  salinity  changes  and  was  calculated  to  be 
2.01,  and  the  resultant  theoretical  flux  ratio  using  the 
equation  of  Radko  and  Smith  is  0.53.  The  line  of  this  slope 
is  shown  in  Figure  23  below  each  set  of  layer  points  to 
show  the  correlation  between  the  flux  ratio  and  lateral 
density  ratio. 
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Figure  23.  Temperature-salinity  values  taken  from  three 
individual  layers  in  the  model .  Each  color 
represents  a  separate  layer. 

The  lateral  coherence  is  clearly  displayed  in  Figures 
24  and  25,  which  show  a  series  of  vertical  temperature 
profiles  from  the  numerical  model  taken  in  the  C-SALT  area 
at  horizontal  separation  of  10  km  in  the  east-west 
direction,  covering  a  total  distance  of  290  km.  The 
temperature  scale  is  correct  for  the  far  left  profile,  with 
each  consecutive  profile  offset  by  1.6  °C  .  Individual 
layers  can  be  traced  along  the  entire  distance,  displaying 
a  very  clear  lateral  coherence,  which  extends  over  300  km. 
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Figure  24.  Series  of  offset  temperature  profiles  from  the  DNS 
in  an  east  to  west  section  with  10  km  separation.  The 
right-most  profile  lies  to  the  east  of  the  others; 
the  total  distance  covered  is  290km.  The  temperature 
scale  is  correct  for  the  left-most  profile  and  each 

subsequent  profile  is  offset  1.6  C. 


Figure  25.  Temperature  profiles  from  Figure  24,  focused  at 

the  depth  of  staircase  formation. 

33 


4 .  Alignment  of  Interface  Slopes  with  Background 

Flow 

To  attempt  to  correlate  the  slope  of  the  interface 

with  the  background  flow  as  discussed  in  Chapter  II,  the 
background  flow  must  first  be  found.  A  Matlab  code  was  run 
to  identify  the  lateral  temperature  and  salinity,  Tt  and  St , 
along  the  direction  of  the  flow  using  the  following 
formulation,  where  u  and  v  are  velocity  in  x  and  y, 

respectively,  in  a  Cartesian  grid: 

u  v 

T,  =  T=  +  Tv]^  (9) 

|v|  |v| 

(10) 

The  slope  of  the  isotherm  and  isohalines  can  then  be 

defined  as: 

Ti 

tg(aT)  =  — 

T 

Z 

tg(as)  =  f- 

(ii) 

The  mean  background  flow  was  found  to  be  predominantly 
West  South  West  in  the  DNS  throughout  the  C-SALT  area  of 
interest.  The  intent  was  to  determine  the  slope  of  the 
interface  by  identifying  the  layer,  and  then  determining 
the  gradient  change.  Unfortunately,  as  can  be  viewed  in  the 
lateral  coherence  figures,  the  depth  change  of  the  layers 
was  exceptionally  small  over  several  hundred  kilometers 
(layers  are  almost  completely  horizontal) .  This  results  in 
a  zero  degree  angle  for  the  interface,  therefore  the 
following  explanation  of  the  continued  analysis  is  provided 
for  future  research  purposes. 
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Using  Wall's  formulation  for  the  interface  slope 
(Equation  7),  the  theoretical  value  can  be  compared  with 
the  numerical  simulation.  Taking  Equations  9  and  10  and 
shifting  90  degrees  to  the  right  to  evaluate  normal  flow 
gives : 


v  (-u) 

T  =r  +  T 


V 


v 

S  =r  +  S 


vl 

(-u) 


V 


V 


(12) 


where  the  numerator  in  Equations  13  and  14  is  now  given  to 
be  Tn  and  Sn  ,  respectively.  It  is  proposed  that  this  can 
also  be  compared  with  the  interface  angle  to  demonstrate  a 
lack  of  correlation  between  interface  slope  and  isotherms 
or  isohalines  normal  to  the  direction  of  flow. 
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IV.  STAIRCASE  EFFECTS  ON  ACOUSTIC  PROPOGATION 


A .  BACKGROUND 

The  effect  of  thermohaline  staircases  has  previously 
been  studied,  both  within  the  C-SALT  area  and  in  other 
areas  of  the  ocean.  Chin-Bing  et  al .  (1994)  used  data  that 
was  collected  from  the  northeast  coast  of  South  America 
adjacent  to  the  C-SALT  area  in  1985.  These  staircases  were 
formed  as  a  result  of  salt-fingering.  The  study  focused  on 
the  effect  of  the  sound  velocity  profile  as  a  result  of  the 
staircases,  and  resultant  acoustic  propagation  simulations. 
The  study  only  investigated  the  effect  at  50  Hz,  which 
creates  an  acoustic  wavelength  of  similar  size  to  the  step 
size.  A  conclusion  from  the  study  was  that  for  significant 
differences  in  prediction  to  be  observed  both  the  source 
and  receiver  needs  to  be  placed  within  the  depth  range  of 
the  staircases. 

Wall  (2007)  built  on  this  investigation  of  acoustic 
propagation  using  a  parabolic  equation  (PE)  model  for  a 
stepped  profile  taken  from  C-SALT.  The  PE  model  was  run 
using  two  different  source  depths  (one  within  and  one  below 
the  staircase  region) ,  and  three  different  source 
frequencies.  Contrary  to  Chin-Bing' s  findings.  Wall  found 
that  the  source  and  receiver  did  not  need  to  be  at  the  same 
depth  as  the  stepped  layer  for  significant  differences  to 
be  observed  between  a  smooth  and  stepped  profile  acoustic 
prediction.  This  finding  was  similar  to  one  found  in  the 
past  by  Wilson  (2007),  who  studied  propagation  through  a 
stepped  structure  resultant  from  diffusive  convection.  Wall 
also  found  that  as  the  frequency  of  the  source  is  increased 
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the  effect  of  the  staircase  on  the  acoustic  prediction  is 
amplified.  A  limitation  to  these  observations  was 
frequency,  where  1200Hz  was  the  highest  frequency  studied. 

This  study  continues  to  build  on  Wall's  observation, 
using  a  different  model  (Bellhop)  to  compare  past  results, 
and  increasing  the  frequency  to  determine  if  the  decrease 
in  wavelength  continues  to  have  an  impact  on  difference  in 
acoustic  prediction  between  smooth  and  stepped  profiles 
taken  from  C-SALT.  The  analysis  of  the  source  within  the 
staircase  is  investigated  further  than  the  previously 
mentioned  studies,  such  that  the  source  is  set  at  a  depth 
corresponding  to  a  prominent  interface,  and  also  at  the 
center  of  a  prominent  mixed  layer  to  determine  whether  the 
positioning  of  the  source  within  the  staircase  region  plays 
a  pivotal  role  in  the  propagation.  Based  on  the  findings  of 
lateral  staircase  alignment  in  earlier  chapters,  an 
acoustic  prediction  is  carried  out  over  a  range  dependent 
simulation  to  assess  how  laterally  coherent  staircases 
impact  acoustic  propagation. 

B .  PROCESS 

The  Bellhop  model  is  a  beam  tracing  model  for  ocean 
acoustic  predictions.  The  theory  behind  the  ray  equation 
used  in  Bellhop  can  be  expressed  as: 

dr  „ 

—  =  c£(s) 
ds 

d£  1  dc 

=  (13) 

V  =  c^(s) 

ds 

d^  _  1  dc 

ds  c2  dz 
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where 


r  ( s )  and  z  ( s ) 


represent  ray  coordinates 


m 


cylindrical  coordinates,  s  is  the  arc  length  along  the  ray, 

and  c(s)  [  ^(s) ,  £(s)  ]  represents  the  tangent  versor  along  the 

ray.  Initial  conditions  for  r(s),  z  (s)  ,  £(s)  and  ^(s)  are 

r  (0)  =r 
z  (0)  =  zs 

e 

_ s_ 

(14) 

„  sin  0 

C(  0)  = - ^ 

c 

s 

where  0s  represents  launching  angle,  (r, zs )  is  the  source 
position,  and  cs  is  sound  speed  at  source  position.  The 
coordinates  are  sufficient  to  obtain  the  ray  travel  time: 


5(0)  = 


cos 


C  ds 
|  c(s) 


(15) 


calculated  along  the  curve  [r(s),  z (s) ]  (Rodriguez,  2012). 

To  run  the  desired  profiles  through  the  Bellhop  model 
an  environmental  file  first  needs  to  be  created.  Data  from 
a  C-SALT  CTD  (Conductivity  Temperature  Depth)  probe  from 
station  20  ( 14°12.70  N  56°17.60  W)  that  displayed  significant 

staircase  development  was  used,  as  was  the  MITgcm  modeled 
simulation  of  laterally  coherent  well-developed  staircases, 
the  temperature  profile  of  these  examples  is  in  Figure  26. 
The  sound  speed  profile  (SSP)  was  calculated  using  MATLAB 
code  from  the  seawater  toolbox,  and  is  displayed  in  Figure 
27.  The  small  increase  in  sound  speed  at  the  base  of 
individual  layers  in  the  staircase  region  SSP  is  expected 
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due  to  the  isothermal  feature,  and  increasing  depth  (and 


hence  pressure),  of  the  layers  within  the  staircase. 
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Figure  26.  C-SALT  Station  20  Temperature  profile  (left)  with 

upper  1000m  displayed  (right) . 
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Sound  Velocity  (m/s) 


Figure  27.  C-SALT  Station  20  SSP  (left)  with  upper  1000  m 

displayed  (right) . 


The  C-SALT  SSP  was  also  smoothed  to  remove  the  sharp 
changes  in  sound  speed  due  to  the  staircase,  using  a  low 
pass  moving  average  filter  with  a  span  of  41  points,  and 
this  profile  was  used  as  the  example  of  an  ocean 
environment  without  staircases,  and  is  in  Figure  28  and  29. 
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Sound  Velocity  (m/s) 


28.  C-SALT  Station  20 


Sound  Velocity  (m/s) 

Figure  29.  C-SALT  Station  20 

smoothed  SSP 


SSP  (left)  and  smoothed  profile 
(right) . 


SSP  to  depth  700  m  (left)  and 
to  depth  700  m  (right) 
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When  building  the  environment  file,  only  1000  m  of 
data  was  used,  with  the  bottom  type  acting  as  an  acoustic 

half-space,  with  a  1.8  g /cm  density  and  attenuation 
coefficient  of  2  dB/m.  The  attenuation  coefficient  was 
chosen  using  Figure  17  of  Hamilton  (1980),  to  be  just  above 
the  attenuation  value  expected  for  compressional  waves 
versus  frequency  in  natural,  saturated  sediments  and 
sedimentary  strata,  for  the  frequency  range  which  we  are 
focused  on  (up  to  5  kHz)  .  The  reasoning  behind  this  was  to 
absorb  as  much  sound  as  possible  into  the  "false"  bottom, 
so  as  to  reduce  the  impact  reflected  or  refracted  rays  re¬ 
entering  the  staircase  region  from  the  bottom  may  have  on 
the  acoustic  propagation. 

The  scenarios  presented  in  part  b  are  conducted  over 
three  source  depths:  5  m,  387  m  and  450  m.  5  m  is  a  typical 
depth  of  a  fixed  hull  mounted  sonar,  and  thus  provides  a 
good  representation  of  a  source  outside  of  staircase 
depths.  378  m  is  the  depth  at  which  a  rather  prominent 
layer  interface  occurs  in  the  C-SALT  profile,  and  475  m  is 
close  to  the  vertical  center  of  a  layer  within  the 
staircase  region.  All  source  depth  scenarios  were  evaluated 
at  frequencies  of  50  Hz,  2000  Hz,  and  5000  Hz  to  evaluate 
the  effect  of  the  deep  layers  on  smaller  wavelengths. 

C .  RESULTS 

1 .  Source  Depth  5  m 
a .  Ray  Trace 

With  both  the  smoothed  and  stepped  profile  the 
ray  trace  (Figure  30)  appears  the  same,  regardless  of 
frequency.  This  is  largely  to  be  expected  as  both  smoothed 
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and  stepped  scenarios  display  the  same  profile  to  a  depth 
of  250  m,  at  which  point  the  staircase  is  introduced.  Thus 
a  source  near  the  surface  experiences  surface  ducting  due 
to  the  mixed  layer,  and  then  a  downward  refracting  ray 
trace  due  to  the  decreasing  sound  speed  profile,  with  a 
shadow  zone  between  the  duct  and  the  downward  refracted 
rays.  The  staircase  has  little  impact  on  this  ray  trace,  as 
the  layer  depth  is  not  large  enough  to  result  in  a 
significant  refractive  change  in  the  overall  downward  ray 
trace  pattern. 


Range  (m) 

Figure  30.  Ray  Trace  profile  for  a  source  at  5  m  in  a 


staircase  environment.  The  ray  trace  is  nearly 
identical  to  that  which  occurs  for  a  smoothed 

profile . 

Jb.  Transmission  Loss 

At  a  source  depth  of  5  m,  there  was  no  difference 
in  the  transmission  loss  between  a  stepped  and  smoothed 
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profile.  The  results  were  the  same  as  frequency  increased. 
Figures  31  and  32  show  the  transmission  loss  experienced  by 
both  profiles  at  2000  Hz.  This  result  is  contrary  to  the 
findings  by  both  Wilson  (2007)  and  Wall  (2007),  which 
stated  the  source  could  be  outside  of  the  staircase  region 
and  have  a  significant  effect  on  propagation. 


Range  (m) 


Figure  31.  Coherent  transmission  loss  (dB/m)  for  a  source  at 
5  m  depth  and  2000  Hz  frequency  for  a  smoothed 

environment . 
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Figure  32.  Coherent  transmission  loss  (dB/m)  for  a  source  at 
5  m  depth  and  2000  Hz  frequency  for  a  stepped 

environment . 


2 .  Source  Depth  387  m 
a .  Ray  Trace 

The  ray  trace  for  a  source  at  387  m  is 
significantly  different  between  the  smooth  and  stepped 
profile  as  seen  in  Figures  33  and  34.  At  this  depth,  the 
source  lies  at  an  interface  for  the  stepped  profile,  but  in 
the  smoothed  profile  at  this  depth  the  SSP  continues  to 
decrease.  It  is  evident  that  when  the  source  lies  at  the 
same  depth  as  an  interface  in  a  region  of  staircases,  there 
is  a  significant  divergence  of  the  rays  at  this  point.  For 
low  launch  angles  the  rays  refract  upward  or  downward, 
leaving  a  rather  significant  shadow  zone  horizontally  in 
front  of  the  source. 
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Figure  33. 


Ray  Trace  for  source  at  387  m  (depth  of  an 
interface)  in  a  staircase  environment. 


Figure  34 


0  1000  2000  3000  4000  5000  6000  7000  8000  9000  10000  11000 
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Ray  trace  for  a  source  at  387  m  in  a  smoothed 

profile . 
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Jb.  Transmission  Loss 

Figures  35  to  40  displays  transmission  loss  for 
both  stepped  and  smooth  profiles  when  the  source  is 
positioned  at  an  interface,  from  50  to  5000  Hz.  As  was 
evidenced  in  the  ray  trace  of  the  stepped  profile,  a  large 
shadow  zone  exists  at  this  depth  at  all  frequencies  tested, 
which  is  not  observed  in  the  smooth  profile  (although  at 
some  frequencies  there  appears  to  be  a  small  area  of 
increased  transmission  loss  above  this  in  the  smoothed 
profile)  .  There  is  also  an  apparent  weak  concentration  of 
signal  in  the  area  just  below  this  shadow  zone,  reasons  for 
this  are  unexplained. 


Range  (m) 

Figure  35.  Coherent  TL  (dB/m)  for  stepped  profile  with 

source  in  interface  at  50  Hz 
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Figure  36. 
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Coherent  TL  (dB/m)  for  smooth  profile  with  50  Hz 

source  at  387  m. 


Range  (m) 


Figure  37.  Coherent  TL  (dB/m)  for  stepped  profile  with 

source  in  interface  at  2000  Hz. 
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Figure  38 
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Coherent  TL  (dB/m)  for  smooth  profile  with  2000 
Hz  source  at  387  m. 


Increasing  the  frequency  to  2000  Hz  appears  to  have  a 
pronounced  effect  on  the  shadow  zone  in  the  stepped 
profile . 
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Figure  39.  Coherent  TL  (dB/m)  for  stepped  environment  with 

source  in  interface  at  5000  Hz. 
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Figure  40.  Coherent  TL  (dB/m)  for  5000  Hz  source  in  smooth 

profile  at  387  m. 

There  is  little  obvious  difference  between  2000  Hz  and 
5000  Hz  when  the  source  is  placed  at  the  interface  of  a 
staircase 

3.  Source  Depth  450  m 
a .  Ray  Trace 

In  this  scenario  the  source  is  positioned  at  a 
depth  that  corresponds  to  the  center  of  a  vertically  mixed 
layer  in  the  stepped  profile.  Unlike  the  previous  example, 
in  this  case  the  low  angle  rays  converge,  causing  a  very 
obvious  duct  to  occur  at  the  source  depth.  In  the  smoothed 
profile,  the  ducting  is  not  evident  to  the  same  extent, 
with  no  particular  features  of  convergence  or  divergence 
appearing  to  stand  out. 
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Depth  (m)  [3.  Depth  (m) 


Range  (m) 


gure  41.  Ray  trace  with  source  at  450  m  (center  of 

vertically  mixed  layer)  in  a  staircase  environment. 


Range  (m) 

Figure  42.  Ray  trace  with  source  depth  450  m  in  a  smoothed 

profile . 
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Jb. 


Transmission  Loss 


The  ducting  that  is  observed  in  the  ray  trace 
using  the  stepped  profile  for  a  source  at  the  center  of  the 
layer  is  also  observed  in  all  the  TL  plots.  In  the  smoothed 
profile,  all  the  plots  are  void  of  any  ducting 
characteristics  resultant  from  staircases.  This 
demonstrates  that  if  the  source  is  placed  within  a  layer  of 
the  staircase,  that  much  higher  horizontal  acoustic  ranges 
may  be  achieved  than  if  there  were  no  staircase  present. 
Once  again,  examples  are  shown  with  lower  and  higher  source 
frequencies.  It  is  evident  that  an  increase  in  frequency 
from  50  Hz  to  2000  Hz  results  in  higher  propagation  ranges 
within  the  duct,  as  a  result  of  increased  leakage  at  lower 
frequencies,  as  opposed  to  trapping  of  smaller  wavelengths. 
However,  there  is  not  a  significant  difference  between  the 
ducted  propagation  of  2000  and  5000  Hz,  and  it  is  proposed 
that  this  is  a  result  of  the  increased  absorption  loss  at 
higher  frequencies. 
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Figure  43.  Coherent  TL  (dB/m)  for  stepped  profile  with  50  Hz 

source  in  layer. 
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Figure  44.  Coherent  TL  (dB/m)  for  50Hz  source  at  450m  in  a 

smoothed  environment. 


54 


0  2000  4000  6000  8000  10000 

Range  (m) 


Figure  45.  Coherent  TL  (dB/m)  for  stepped  profile  with  2000 

Hz  source  in  layer. 
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Figure  46.  Coherent  TL  (dB/m)  for  2000  Hz  source  at  450m  in 

smoothed  profile. 
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Figure  47.  Coherent  TL  (dB/m)  for  5000  Hz  source  in  layer  of 

stepped  profile. 
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Figure  48.  Coherent  TL  (dB/m)  for  5000  Hz  source  at  450  m  in 

smoothed  profile. 
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4. 


Range-Dependent  Source  in  Layer 


a .  Ray  Trace 

Following  on  from  the  investigation  of  a  source 
placed  in  a  staircase  layer  as  opposed  to  a  staircase 
interface,  we  investigate  a  range-dependent  example  of 
acoustic  propagation  in  staircases.  Figure  24  displayed  a 
series  of  temperature  profiles  from  the  numerical  model 
spaced  every  10  km  with  lateral  coherence  over  a  range  of 
300  km.  Nine  adjacent  profiles  from  this  model  were 
converted  to  SSP  using  the  MATLAB  acoustic  toolbox  and  then 
compiled  into  a  range-dependent  input  file  compatible  with 
Bellhop.  Figure  50  shows  the  individual  sound  speed 

profiles,  and  Figure  51  displays  an  interpolated  sound 
speed  over  the  entire  range  showing  the  subtle  shift  in 
sound  speed  over  range  and  depth  as  a  result  of  small 
vertical  shifts  in  the  laterally  coherent  layers.  The 

center  of  a  35  m  vertical  layer  at  a  range  of  0  m  within 
the  staircase  region  was  identified  at  a  depth  of  785  m, 
and  a  ray  trace  was  conducted  with  the  source  at  this 
depth.  Ducting  is  evident  at  this  depth  as  can  be  seen  by 
the  red  ray  path  in  Figure  49.  An  interesting  feature  is  a 
shift  in  the  ducted  rays  at  a  range  of  approximately  60  km 

from  the  source.  It  is  hypothesized  that  this  is  a  result 

of  a  subtle  shift  in  the  depth  of  interfaces  of  connecting 
layers  as  range  is  increased  away  from  the  source.  It  is 
apparent  from  the  ray  trace  that  the  duct  shifts  up  or  down 
to  remain  within  the  layer,  as  the  layer  thickness  changes 
(interfaces  shift  up/down),  over  extended  ranges. 
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Figure  49.  Ray  trace  for  source  at  785  m  (source  in  layer), 
in  a  range  dependent  staircase  environment. 


Figure  50.  Range-dependent 


ound  speed  profiles  at  10  km 
spacing . 
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Figure  51.  Sound  speed  (m/s)  as  a  function  of  depth  and 
range  due  to  laterally  coherent  staircases. 

Jb.  Transmission  Loss 

Figure  52  to  55  displays  transmission  loss  for 
the  range-dependent  case  when  the  source  is  placed  at  785  m 
within  a  prominent  layer,  over  frequencies  of  50  to  5000 
Hz.  A  notable  difference  in  analysis  to  earlier  TL  plots  is 
that  the  TL  is  now  being  observed  over  a  range  of  100  km. 
Ducting  is  evident  at  the  source  depth  in  all  cases  with 
the  exception  of  50  Hz.  In  the  50  Hz  scenario  there  is 
increased  energy  at  this  depth  at  short  range,  as  a  result 
of  leakage.  However,  the  larger  wavelength  means  sound  is 
not  trapped  in  this  layer.  The  best  example  of  ducting 
occurs  at  500  Hz,  where  there  is  very  obvious  trapping  of 
the  rays  through  the  entire  range,  with  significant  energy 
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still  remaining  at  80  km  from  the  source.  The  ducting  is 
also  obvious  at  2000  Hz,  but  is  largely  diminished  by 
5000  Hz.  It  can  be  surmised  that  as  the  frequency  of  the 
source  is  increased  the  absorption  at  the  higher 
frequencies  quickly  becomes  a  more  dominant  process. 
Therefore,  in  a  range-dependent  situation  a  low  to  medium 
frequency  source  of  approximately  1000  Hz  appears  most 
effective  in  propagation  of  sound  within  a  duct.  It  is 
proposed  that  this  frequency  is  highly  dependent  on  the 
thickness  of  the  vertical  layers  also,  as  small  layers  are 
likely  to  only  trap  higher  frequency  rays. 
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Figure  52.  Coherent  TL  (dB/m)  for  50  Hz  source  within  layer. 
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Figure  53. 
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Figure  54.  Coherent  TL  (dB/m)  for  2000  Hz  within  layer. 
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Figure  55.  Coherent  TL  (dB/m)  for  5000  Hz  within  a  layer. 

5 .  Range-Dependent  Source  in  Interface 
a .  Ray  Trace 

Using  the  same  range-dependent  profiles  from  part 
4,  a  source  was  placed  at  a  depth  of  750  m,  which  coincides 
with  an  interface.  In  the  resultant  ray  trace  from  this 
there  is  no  apparent  ducting  of  the  rays.  Similarly  to  what 
was  seen  in  the  range-independent  scenario  previously 
discussed,  there  appears  to  be  a  divergence  of  the  rays  in 
a  horizontal  direction  from  the  sound  source  as  highlighted 
in  the  circle.  Close  to  the  sound  source  this  may  generate 
a  shadow  zone,  however  in  the  range  dependent  case  this 
area  is  quickly  flooded  with  energy  from  refracted  and 
reflected  rays. 
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Figure  56.  Ray  trace  for  source  at  staircase  interface  for 

range-dependent  profiles. 

Jb.  Transmission  Loss 

Transmission  loss  was  again  calculated  from  50  to 
5000  Hz  for  the  source  at  the  interface.  It  appears  that 
some  energy  still  manages  to  get  trapped  within  the  lower 
layer  ducts,  and  there  is  faint  evidence  of  this  at  a  large 
range  from  the  source  of  approximately  70  km,  at  500  and 
2000  Hz  frequencies.  This  ducting  effect  is  much  less 
pronounced  when  the  source  is  at  the  interface  as  opposed 
to  the  center  of  the  layer. 
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Figure  58.  Coherent  TL  (dB/m)  for  500  Hz  source  at 

interface . 
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Figure  60.  Coherent  TL  (dB/m)  for  5000  Hz  source  at 

interface  depth. 
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V. 


CONCLUSIONS 


Thermohaline  staircases  are  a  significant  feature  of 
ocean  mixing,  and  this  study  has  aimed  to  explore  staircase 
evolution  in  three  dimensions  as  a  result  of  salt 
fingering,  through  numerical  simulation. 

A  large  task  was  undertaken  in  attempting  to  simulate 
three  dimensional  staircases  over  an  entire  ocean  basin.  To 
the  authors'’  best  knowledge  this  is  the  first  basin-scale 
simulation  that  exhibits  spontaneous  layering.  A  promising 
result  has  been  the  modeling  of  laterally  coherent 
staircases  over  several  hundreds  of  kilometers  in  the 
ocean.  This  provides  an  interesting  area  of  further 
research,  particularly  with  the  alignment  of  layer 
interfaces  with  background  isotherms.  While  this  was  not 
achieved  in  this  thesis,  it  is  proposed  that  Wall's 
hypothesis  of  a  correlated  layer  slope  with  thermohaline 
will  hold  true,  and  with  more  computing  hours  and  time,  it 
is  believed  that  a  high  resolution  grid  would  make  this  an 
achievable  task  for  researchers  in  the  future. 

This  study  aimed  to  parameterize  turbulent  diffusivity 
through  the  process  of  trial  and  elimination  in  order  to 
closely  simulate  the  thermohaline  staircases  seen  in  the  C- 
SALT  observation  area.  Through  this  process,  the  turbulent 

diffusivity  was  indirectly  estimated  at  2.5  x  10~6  m2  /  s  . 

A  comparison  between  the  gamma-effect  variation  in  the 
flux  ratio  with  the  density  ratio  confirmed  a  very  distinct 
link  between  the  flux  ratio,  the  density  ratio,  and 
staircase  formation.  This  lends  further  support  to  the  flux 
ratio  model  of  Radko  and  Smith  (2012) . 
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The  numerical  simulation  of  merging  events  indicates 
the  prevalence  of  instability  that  leads  to  the  'B  type 
merger'  as  discussed  in  Radko  (2007),  and  this  simulation 
clearly  supports  the  theory  that  development  of 
thermohaline  staircases  by  salt  fingers  is  a  direct  result 
of  strong  interfaces  growing  at  the  expense  of  weaker  ones 
which  ultimately  erode. 

The  acoustic  propagation  modeling  through  the 
staircase  of  both  the  C-SALT  profiles  and  the  numerical 
model  indicate  that  thermohaline  staircases  have  a  very 
pronounced  effect  on  sound  propagation,  which  is 
accentuated  with  position  placement  of  the  sound  source 
within  the  staircase.  The  ducting  that  occurs  when  a  source 
is  placed  centrally  within  a  layer  of  a  staircase  is  a 
significant  feature,  and  the  increased  range  of  the  duct 
with  lateral  coherence  of  staircases  is  an  extremely 
important  facet  to  be  considered  in  acoustic  modeling. 
Understanding  where  staircase  regions  of  lateral  coherence 
occur,  coupled  with  knowledge  of  interface  and  layer  depths 
within  this  region  can  be  of  critical  importance  for 
tactical  employment  of  submarines,  gliders,  variable  depth 
sonars,  torpedo  homing,  and  other  underwater  vehicles 
within  this  region.  This  information,  combined  with  an 
awareness  of  the  impact  of  differing  source  frequencies  on 
the  ducting,  can  result  in  a  significant  exploitation  of 
the  environment  for  tactical  purposes. 
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